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ABSTRACT 

Models for the evolution of the integrated energy spectrum of primary cosmic ray electrons in clusters of 

■ galaxies have been calculated, including the effects of losses due to inverse Compton (IC), synchrotron, and 
| bremsstrahlung emission, and Coulomb losses to the intracluster medium (ICM). The combined time scale for 

t-H i these losses reaches a maximum of ^3 x 10 9 yr for electrons with a Lorentz factor 7 ~ 300. A variety of models 

■ for the time evolution of particle injection are considered, including models in which the electrons are all produced 
at a single epoch in the past, models with continuous particle acceleration, and combinations of these. Analytical 

|"T ^ • solutions are given for a number of limiting cases. Numerical solutions are given for more general cases. Only 

clusters in which there has been a substantial injection of relativistic electrons since z <J 1 will have any significant 
population of primary cosmic ray electrons at present. For models in which all of the electrons were injected in 
the past, there is a high energy cutoff -f max to the present electron distribution. At low energies where Coulomb 
losses dominate, the electron distribution function N(j) tends to a constant value, independent of 7. On the other 
hand, if electrons are being accelerated at present, the energy distribution at high and low energies approaches 
' steady state. If the electrons are injected with a power-law distribution, the steady state distribution is one power 

\Q [ steeper at high energies and one power flatter at low energies. In models with a large initial population of particles, 

but also with a significant rate of current particle injection, the electron distributions are a simple combination the 
t-H ' behavior of the initial population models and the steady injection models. There is a steep drop in the electron 

population at 7 mflv , but higher energy electrons are present at a rate determined by the current rate of particle 
injection. Increasing the ICM thermal gas density decreases the number of low energy electrons (7 < 100). If the 
magnetic field is greater than generally expected, (B p> 3 /iG), synchrotron losses will reduce the number of high 
energy electrons. A significant population of electrons with 7 ~ 300, associated with the peak in the particle loss 
time, is a generic feature of the models, as long as there has been significant particle injection since z <J L 
1. \ The IC and synchrotron emission from these models was calculated. In models with steady particle injection 

with a power-law exponent p, the IC spectra relax into a steady-state form. At low energies, the spectrum is a 
power-law with a « -0.15, while at high energies a -1 . 15. These two power-laws meet at a knee at v ~ 3 x 10 16 
Hz. In models with no current particle injection, the cutoff in the electron distribution at high energies (7 > 7 mflA ) 
results in a rapid drop in the IC spectrum at high frequencies. In models in which the current rate of particle 
injection provides a small but significant fraction of the total electron energy, the spectra show an extended hump 
at low frequencies (v <; 10 17 Hz), with a rapid fall off above v ~ 10 16 Hz. However, they also have an extended 
hard tail of emission at high frequencies, which has a power-law spectrum with a spectral index of a -1.15. 

In the models, EUV and soft X-ray emission are nearly ubiquitous. This emission is produced by electrons 
with 7 ~ 300, which have the longest loss times. The spectra are predicted to drop rapidly in going from the 
EUV to the X-ray band. The IC emission also extends down the UV, optical, and IR bands, with a fairly flat 
spectrum (-0.6 ;> a ;> +0.3). At hard X-ray energies ;>20 keV, IC emission should become observable against 
the background of thermal X-ray emission. Such hard X-ray (HXR) emission is due to high energy electrons 
(7 ~ 10 4 ). The same electrons will produce diffuse radio emission (cluster radio halos) via synchrotron emission. 
Because of the short loss times of these particles, HXR and diffuse radio emission are only expected in clusters 
which have current (or very recent) particle acceleration. Assuming that the electrons are accelerated in ICM 
shocks, one would only expect diffuse HXR/radio emission in clusters which are currently undergoing a large 
merger. The luminosity of HXR emission is primarily determined by the current rate of particle acceleration in 
the cluster. The spectra in most models with significant HXR or radio emission are approximately power-laws 
with a w — 1.1. The IC spectrum from the EUV to HXR are not generally fit by a single power-law. Instead, there 
is a rapid fall-off from EUV to X-ray energies, with a power-law tail extending into the HXR band. 

Subject headings: cosmic rays — galaxies: clusters: general — intergalactic medium — radiation mechanisms: 
nonthermal — ultraviolet: general — X-rays: general 

1. INTRODUCTION to hot (~ 10 8 K), diffuse intracluster gas (see Sarazin [1988] 

r . v , .. jj, . . , , for a review). This intracluster medium (ICM) contains heavy 

Early X-ray observations with the uhuru satellite established , • , , , , • , , r , 

. , . c . . . . . a a V elements with abundances, relative to hydrogen, of about one 

that clusters of galaxies were luminous, extended X-ray sources , . , „ , , , ' • , , , 

:n , ,. o rr , i/v7i\ Tu v * a third of the solar values. The gas is thought to have originated 

(Cavahere, Gursky, & Tucker 1971). The X-ray emission is due . , . . . „ „. , fe . , fe , , , , 

mainly from infall of intergalactic gas, although the heavy ele- 
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ments probably came from gas ejected by galaxies. The gas is 
believed to have been shock-heated during infall into the clus- 
ter, during subcluster mergers, and as a portion of it was ejected 
by galaxies. 

It is useful to compare the intracluster medium with the inter- 
stellar medium which fills the volumes of space between stars 
in our Galaxy. The interstellar medium in our Galaxy contains 
a mixture of at least three different types of matter: thermal 
gases of various densities and temperatures; solid dust grants; 
and relativistic materials, including magnetic fields and cosmic 
ray particles. To what extent would one expect to find similar 
materials in the intracluster medium? It may be difficult for 
dust grains and cooler thermal gas phases to survive in the hot 
ICM. However, there is no obvious reason why clusters should 
not contain magnetic fields and relativistic particles. 

In fact, there are a number of reasons to think that cosmic 
ray particles might be particularly abundant within the intra- 
cluster medium. First, clusters of galaxies should be very ef- 
fective traps for cosmic ray ions and electrons. Under reason- 
able assumptions for the diffusion coefficient, particles with 
energies of less than <;10 6 GeV have diffusion times which 
are longer than the Hubble time (Berezinsky, Blasi, & Ptuskin 
1997; Colafrancesco & Blasi 1998). Second, the lifetimes of 
cosmic ray particles, even the electrons which are responsible 
for most of the radiative signatures of relativistic particles, can 
be quite long. The radiation fields (optical/IR and X-ray) and 
magnetic fields (B <; 1 /iG) in the ICM are low enough that high 
energy electrons mainly lose energy by inverse Compton (IC) 
scattering of Cosmic Microwave Background (CMB) photons 
(Sarazin & Lieu 1998). Lower energy electrons can lose en- 
ergy by Coulomb interactions with the plasma; however, at the 
very low densities (n e <J 10~ 3 cm" 3 ) in the bulk of the ICM, this 
is only important for electrons with 7 < 200 (Sarazin & Lieu 
1998). (Here, 7 is the Lorentz factor for the electrons, so that 
their total energy is jninC 2 .) For electrons with 7 ;> 200, the 
lifetime is set by IC losses and is 
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The lifetimes of ions are set by interactions; for protons, this 
gives tio„ ;> 10 n ("e/10~ 3 cm" 3 ) _1 yr. Thus, clusters of galaxies 
can retain low energy electrons (7 ~ 300) and nearly all cosmic 
ray ions for a significant fraction of a Hubble time. 

Cluster of galaxies are likely to have substantial sources of 
cosmic rays. Clusters often contain powerful radio galaxies, 
which may produce and distribute cosmic rays throughout the 
cluster, in addition to possible contributions from the cluster 
galaxies. However, the primary reason why the cosmic ray pop- 
ulations in clusters might be large is connected with the high 
temperature of the intracluster gas. This indicates that all of 
the intracluster medium (typically, 1O 14 M of gas) has passed 
through strong shocks with shock velocities of ~ 1 000 km/s dur- 
ing its history. In our own Galaxy, whenever diffuse gas under- 
goes a strong shock at velocities of this order, a portion of the 
shock energy goes into the acceleration of relativistic particles 
(Blandford & Ostriker 1978; Bell 1978a,b; Blandford & Eich- 
ler 1987; Jones & Ellison 1991). Thus, it seems likely that 
relatively efficient particle acceleration also occurs in clusters 
of galaxies. 

Direct evidence for the presence of an extensive population 
of relativistic particles and magnetic fields in the ICM comes 
from the observation of diffuse synchrotron radio halos in clus- 



ters (e.g., Wilson 1970; Hanisch 1982; Giovannini et al. 1993; 
Deiss et al. 1997). More recently, extreme ultraviolet (EUV) 
and very soft X-ray emission has been detected from a num- 
ber of clusters (Lieu et al. 1996a,b; Mittaz, Lieu, & Lockman 
1998). Although the origin and even the existence of this radi- 
ation remain controversial, one hypothesis is that it is inverse 
Compton (IC) emission by relativistic electrons (Hwang 1997; 
EnGlin & Biermann 1998; Sarazin & Lieu 1998). Finally, there 
are a number of reports of detections of hard X-ray emission 
from clusters of galaxies with BeppoSAX, which might be due 
to IC emission from higher energy electrons (Fusco-Femiano 
et al. 1998; Kaastra, Bleeker, & Mewe 1998), although this 
interpretation is still uncertain (Henriksen 1998; Goldoni et 
al. 1998). The EUV emission would require electrons with 
7 ~ 300, while the hard X-ray emission would require 7 ~ 10 4 . 

The BeppoSAX detection of hard X-ray emission from Coma 
implies that the average magnetic field over large volumes of 
the cluster is 0.16 ^G. Faraday rotation measurements towards 
Coma and other clusters suggest that the magnetic field in the 
core of the cluster is stronger, perhaps around 1 pG (Kim et al. 
1990; Kim, Kronberg, & Tribble 1991; Kronberg 1994). In our 
models, we will use 0.3 /iG as a typical value for the intracluster 
magnetic field on large scales. 

The energy spectra of relativistic particles in many astro- 
physical environments are often modeled as power-laws. The 
energy spectrum of cosmic ray electrons in clusters of galax- 
ies have been represented as power-laws in most investigations 
(e.g, Hwang 1997; EnGlin & Biermann 1998). This should be a 
reasonable approximation at energies where all of the particles 
are relativistic and where there are no other physical scales af- 
fecting the energy. However, there are a number of reasons why 
the populations of relativistic electrons in clusters of galaxies 
might not be well-represented by power-laws. First, at the low 
energies of interest for the EUV emission (7 ~ 300, or E ~ 150 
MeV, where E is the particle kinetic energy), the electrons are 
relativistic, but the ions are not. Second, at the energies inter- 
est for the EUV emission (7 ~ 300), the life times of the elec- 
trons are close to the Hubble time (eq. [l]). At much higher or 
lower energies, the electron life times are much shorter. Thus, 
losses affect different parts of the electron energy spectrum in 
different ways. Finally, at low energies 7 <J 200, the electrons 
mainly lose energy to the plasma due to Coulomb losses. At 
higher energies 7 > 200, the most important losses are due to 
IC or synchrotron emission. These different loss mechanisms 
have different energy dependences, and this breaks any simple 
scaling in the particle losses. 

In this paper, illustrative models will be calculated for the 
energy spectrum of primary, relativistic electrons. Recently, 
models in which the relativistic electrons are secondaries pro- 
duced by interactions of cosmic ray ions have been given by 
Colafrancesco & Blasi (1998). Cosmic ray ions and secondary 
electrons will be included in a subsequent set of calculations. 
In this paper, we will concentrate on the total energy spectrum 
of all the relativistic electrons in the cluster, under the assump- 
tion that they remain trapped in the cluster (Berezinsky et al. 
1997; Colafrancesco & Blasi 1998). In a subsequent paper, we 
will discuss the spatial distribution of the cosmic ray particles. 
In this paper, radio synchrotron and inverse Compton emission 
spectra will be determined up to the hard X-ray spectral band. 
Gamma-ray emission by the same electrons (as well as emis- 
sion by ions through ir° decay and other nuclear reactions) will 
be included in a later paper. 
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2. EVOLUTION OF THE ELECTRON SPECTRUM 

2.1. Evolution Equation 

The evolution of the cosmic ray electron population in clus- 
ters is given by the diffusion-loss equation (Ginzburg & Sy- 
rovatskii 1964): 

^P-=-n(E)V-y+V-[D(E)Vn(E)]+-^[b(E)n(E)]+q(E), 

(2) 

where n(E)dE is the number density of electrons with kinetic 
energies in the range E to E + dE, v is the velocity of the 
ICM, djdt is the Lagrangian time derivative, D(E) is the dif- 
fusion coefficient, and q(E)dE gives the rate of production of 
new cosmic ray particles per unit volume with energies in the 
range E to E + dE. The losses by the particles are given by 
b(E) = -{dE I dt), where the derivative gives the rate of change 
in the energy of a single particle with an energy of E. This for- 
mulation of the losses assumes that they are continuous; this is a 
reasonable approximation for the Coulomb and inverse Comp- 
ton losses which are most important for electrons in clusters, 
but is not completely correct for bremsstrahlung losses. Equa- 
tion (||) ignores the acceleration of the electrons (which we will 
treat separately and include in q) and their generation as secon- 
daries through the interaction of other particles. 

Here, we consider models for the total cosmic ray electron 
spectrum of clusters, without considering the spatial distribu- 
tion of the particles. Let N(E)dE be the total number of elec- 
trons in the cluster with kinetic energies in the range E to 
E + dE. I will also assume that cosmic ray electrons and ions 
are trapped within clusters, and cannot escape. Models for cos- 
mic ray diffusion within clusters do indeed suggest that they 
should be trapped (Berezinsky et al. 1997; Colafrancesco & 
Blasi 1998). Then, integrating equation (||) over the volume of 
the cluster leads to an equation for the time evolution of the 
total electron spectrum, 



^ = A [KW)]+e(i?) , 



(3) 



where Q(E)dE gives the total rate of production of new cosmic 
ray electrons in the energy range E to E + dE. I have replaced 
b(E) by its value averaged over all the particles in the cluster. 
The first two terms on the right hand side of equation are 
converted by Green's formula into fluxes at the outer surface 
of the cluster, which are zero since we assume that the cluster 
retains all of the particles. 

I will use the Lorentz factor 7 of the electrons as the inde- 
pendent variable rather than the kinetic energy E = (7- l)m e c 2 , 
where m e is the electron mass. With the definitions Nirfd'y = 
N(E)dE, Q(j)dj = Q(E)dE, and 6(7) = b(E)/m e c 2 , the expres- 
sion for the energy loss by an individual particle becomes 



a?7 
dt 



(4) 



and the equation for the evolution of the electron populations is 

(5) 

2.2. Loss Rates 
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The loss function £(7) for inverse Compton scattering of the 
Cosmic Microwave Background (CMB) is given by 



4 (Tj- 2 

bic(l)=T 7 U C mb : 

3 m e c 



1.37 x lo^Va+zyV 1 



(6) 



where oy is the Thomson cross-section, and Uqub is the en- 
ergy density in the CMB at the cluster, and z is the redshift of 
the cluster. The numerical value follows from the present CMB 
temperature of Tqmb = 2.73 K and the redshift evolution of the 
CMB energy density. At very high energies (7 ;> 10 9 ), the IC 
losses are smaller because of lower Klein-Nishina cross-section 
(Gould 1975). 




FIG. 1. — The values of the losses function ^(7) for inverse Compton (IC) 
emission, Coulomb losses, synchrotron losses, and bremsstrahlung losses as a 
function of 7. The values assume n e = 10 -3 cm -3 , B = 1 /iG, and z = 0. 



The expression for the loss rate due to synchrotron radiation 
is the same as that in equation (||) if one substitutes the energy 
density in the magnetic field Ub = B 2 / (871- ) for the CMB energy 
density: 



V>(7) = \— 1 2 Ub = 1.30 x 10" 21 7 2 
3 m e c 



B 
T^G 



(7) 



where B is the magnetic field in the ICM. Both equations (Q) 
and (0) assume 7 3> 1. Obviously, the loss rates for IC and for 
synchrotron losses vary in the same manner with 7. The intra- 
cluster field is estimated to have a values of roughly B ~ 1 /iG 
(e.g., Rephaeli, Gruber, & Rothschild 1987; Kim et al. 1990). 
As long as the field is not much stronger than this, IC losses 
will dominate over synchrotron losses, with 



bsyn _ Ub 
bic Uqmb 



■■ 0.095 (1+z) 4 



1/xG 



(8) 



The relativistic electrons will also lose energy by interactions 
with the thermal plasma. The Coulomb losses due to collisions 
with thermal electrons give a loss rate which is approximately 
(e.g., Rephaeli 1979) 



bcouiil) « 1.2x10- 



1.0 + 



ln(7/n e ) 



75 



(9) 



where n e is the thermal electron density in the cluster. The 
collisions between cosmic ray electrons and thermal ions and 
electrons also produce radiation through bremsstrahlung. For 
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the completely unscreened limit which is appropriate to a low 
density plasma, the loss rate due to bremsstrahlung is given ap- 
proximately by (e.g., Blumenthal & Gould 1970) 



b b re m (D~ 1.51 x 10- 1 V7[ln(7) + 0.36]s" 



(10) 



The numerical coefficient includes electron-ion and electron- 
electron bremsstrahlung, and assumes a helium abundance by 
number of 10% of hydrogen. 

The values of these loss functions as a function of 7 are given 
in Figure [l] for an intracluster medium (ICM) electron density 
of n e = 1.0 x 10~ 3 cm" 3 and a typical ICM magnetic field of 
B = 1 [j,G. It is clear from this Figure and the expressions for 
the loss functions that IC losses are dominant at large energies, 
7 > 200, while Coulomb losses dominate for sufficiently small 
7 <J 200 or for higher densities. Bremsstrahlung losses are un- 
likely to be dominant unless the density is higher than is typical 
in the bulk of the ICM (outside of cooling flow regions). Simi- 
larly, synchrotron losses are unlikely to be dominant unless the 
magnetic field is much stronger than 1 /iG. 

2.3. Electron Lifetimes 
One can define an instantaneous time scale for particle losses 



by 
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Values for this loss time scale at the present time (z = 0) are 
shown in Figure |[ The solid curve gives values assuming an 
average electron density of n e = 10~ 3 cm" 3 and a magnetic field 
of B = 1 /iG. For values of the magnetic field this small or lower, 
synchrotron losses are not very significant, and ti oss is nearly 
independent of B. The short dashed curve shows the effect of 
increasing the magnetic field to B = 5 /iG; the losses at high 
energies are increased, and the loss time scales shortened. The 
dash-dot curve shows the loss time scale if the electron density 
is lowered to n e = 10~ 4 cm" 3 . This reduces the losses at at low 
energies, and increases the loss times there. 



the maximum loss time is quite long, unless either the mag- 
netic field or electron density are larger than typical values for 
most of the volume of a cluster. The maximum loss times are 
~ 3— 10 Gyr, which is comparable to the likely ages of clusters. 
Thus, particl es ca n accumulate with 7 ~ 300 for long periods in 
clusters. In § 3.1 below, we will show that although electrons at 
7 ~ 300 can accumulate for a significant fraction of the Hubble 
time, the increase in the energy density of the CMB implies that 
the accumulation period is restricted to z < 1 . 

3. ANALYTIC SOLUTIONS 

A number of simple analytic solutions for the evolution of 
the cosmic ray population (eq. |5]) are possible. I will consider 
two simple classes of solutions: solutions with an initial popu- 
lations of particles N(j,ti) at some initial time f,- or redshift z; 
but no subsequent injection of additional particles [(2(7) = 0]; 
and solutions with no initial population [iV(7,f,-) = 0] but with 
continual injection of particles at a rate 2(7,0- Since equa- 
tion (Q) is linear in N(j,t), any solution can be written as a 
superposition of these two solutions. 

3.1. Initial Injection Only 

3.1.1. General Solution 

First, consider a case where there is no continual source of 
new particles (Q = 0). Then, the population simply evolves due 
to the loss of energy by individual particles. At this point, as- 
sume that the loss function b(j,t) is an arbitrary function of 7 
and time t . Given an initial value for the energy of a particle 7, 
at time f,-, equation (Q) can be integrated to give the value of 7 
at a later time t. Alternatively, equation (Q) can be integrated 
backward in time to give the initial energy 7,(7, ?,?;) a * U corre- 
sponding to a particle with energy 7 at a later time t. If there 
is no subsequent injection of particles, then the number of par- 
ticles is conserved if one follows their energies. This implies 
that 



r CO rC 

/ N(j',t)dj'= / 



N{i,tddi. 



(12) 



n e = 10 cm 




The differential population density is then given by 



N(j,t) = N(rfi,td 



9ji_ 
c?7 



(13) 



3.1.2. Constant Loss Function 

For periods of time which are short compared to the Hubble 
time, the total loss function £(7) may be approximately inde- 
pendent of time. If a particle starts (at time t = f,) with a Lorentz 
factor 7i, after a time t - f, its energy will have been reduced to 
7, where 

n- di 



■(t-ti). 



(14) 



- The solid curve gives the instantaneous loss time scale ti os , 
(eqn. lilt) as a function of 7 for electrons in a cluster with an electron den- 
sity oin e = 10 -3 cm -3 and a magnetic field of B = 1 ^tG. The short-dash curve 
5 fj,G, while the dash-dot curve is for n e 



sity 
is for B 
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The characteristic features of ti oss for parameters relevant to 
clusters are that it has a maximum at 7 ~ 100-500, and that 



Applying this energy loss to the entire electron population, one 
finds that the electron spectrum evolves as 



N{ 1 ,t) = N[ li ( 1 ,t),t l ] 



17,(7, 01 



(15) 



where 7,(7,0 is given implicitly by equation (|l4j). 
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3.1.3. Solution at High Energies: Non-Cosmological 

IC losses dominate at high energies. If we consider the evo- 
lution over time scales which are shorter than the age of the 
Universe, so that the redshift z is approximately constant, then 
the loss function at high energies can be written as 



7> 10 2 



where the numerical constant b\ is given in equation 
single particle loss equation ( |l4| ) gives 



7> : 



1- 



(16) 
The 

(17) 



After the passage of a time At = (f-f,), all of the electrons with 
energies above 7 m£ « will have been removed, where 



1 



jmax ~ 



biAt 



(18) 



Thus, the resulting electron distribution is related to the initial 
distribution by 



iV( 7 ,0 = 



iV[7/(l-7/7„,„),f,-] 





7 <ln, 
7 > In, 



(19) 



If the initial particle distribution was a power-law with 
N(ji,ti) = N\'yJ Pa , then the resulting final distribution is 



N{ ll t) = 



^ l7 -po (1 _ 7 / 7mflt) Po-2 7<7(? 

7 > 7„ 



(20) 



These equations also apply if there if a significant magnetic 
field in the cluster, since the energy losses from synchrotron 
emission have the same dependence on 7 as those for IC scatter- 
ing. The value of b\ would be increased to include synchrotron 
losses according to equations (||) and (0). 

3.1.4. Solution at High Energies: Cosmological 

These expressions do not apply to the evolution of the elec- 
tron energy spectrum over cosmological time scales because of 
the variation in the energy density of the CMB with redshift z. 
As a result, the loss rate due to IC emission varies with electron 
energy and redshift as (eq. n) 



b IC ( 1 ,z) = b 11 2 (l+z) 4 . 



(21) 



The equation for the evolution of the energy of a single particle 
subject to IC losses is then 



T 2 " 



-bi{\+zfdt. 



(22) 



The age of the Universe t and redshift z are related implicitly 
by the equation of cosmic dynamics (e.g., Weinberg 1972), 



1 

H~o 



(i+zf 



7 _l. 2< ?° 

■2q + 

x 



-1/2 



dx, 



(23) 



here Ho is the Hubble constant and go is the deceleration pa- 
rameter at the present time. 



After some algebra, equation ( |22| ) can be integrated to yield 
the current 7 for an electron which had an initial energy of 7, at 
a redshift zf. 

1 1 bi 

= 77- /(«<>, z,0, (24) 

7 7i H o 

where the function f(qo,z) is defined by 
1 



f(qo,z) = 



I5ql 



For qo = 0, the limit of equation d25|) is 



f(Q,z) = z[l+z+ 



\/l+2q z (2 - 1 0q +I5ql~ IqoZ 
+10^z+3^z 2 )-2+10 ?0 -15^] .(25) 



(26) 




FIG. 3. — The present values of the maximum remaining electron Lorentz 
factor fmax is shown as a function of the redshift z, at which the electrons were 
last injected in the cluster. Only IC losses are included. Curves are shown for 
three cosmological models with values of the cosmic deceleration parameter of 
qo = 0, 0.2, and 0.5. The values of ■y max are scaled to the value of the present 
Hubble constant, with h = (H a /W0 km s~' Mpc~'). 



The cosmological expression for the maximum remaining 
energy of electrons subject to IC losses is then 



]max ~ 



H 



1 



b\ f(qo,zd 



(27) 



Figure || shows values of the maximum Lorentz factor as a 
function of the redshift Zi for cosmological models with val- 
ues of go = 0, 0.2, and 0.5. The values of j max are scaled by 
/^(//o/lOOkms" 1 Mpc" 1 ). 

3.1.5. Solution at High Energies with Synchrotron Losses 

If the ICM magnetic field is much stronger than is currently 
believed, synchrotron losses might be important for high energy 
electrons. Then, the expressions for the cosmological evolution 
of an initial electron population at high energies would be mod- 
ified in a way which would depend on the temporal variation of 
the magnetic field. If the magnetic field were constant, then the 
losses at high energies could be written as 



M7,z)«7 2 [b[ c (l+z) 4 +bf n ] 



(28) 
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where the values of b'f and bf can be found from equations (^) 
and (0). Then the evolution of the energy of a single particle at 
high energies is given by 



1 1 

7 7/ 



ffn 



{b[ C f{q ,zd + bT[g(q Q ,Q)-g(q Q ,Zi)]} , (29) 



where g(qo,z) is the integral on the right side of equation (£3j). 
This function is most easily given in the standard parametric 
form (e.g., Weinberg 1972) 



#(<?(), z) = 



(1+z)" 1 

g0 [sinh(9)-e] 
(1-2 9 „VV 2 



|(i+zr 3/2 



cosh(6>) = 1 + 



l-2go J_ 
90 1+z 



9o = 

0<<7„<l/2 
?o = l/2 



g0 [e-sin(9)] 
(250-1) 3 / 2 



cos(6) = 1 - 



2 g0 -l 1 
qo 1+z 



1/2 < qo 

(30) 

Then, the upper cutoff in the electron energy distribution is 
given by 



//<> 



b[ c f(qo ,zd + bf" [g(q , 0) - g(q ,Zi)l 



(31) 



The resulting particle distributions are still given by equa- 
tions ( |l9| ) or (po|), with the correct value for ~f max . If one wants 
to calculate the particle distribution at an intermediate redshift 
z (0 < z < zi) rather than at present (z = 0), then one substi- 
tutes [/(40vZi)-/^o,zllfor f(<l!hZi) and g(q ,z) for g(q ,0) in 
equations (24), (El), (E9h, and (pi|). 



3.1.6. Solution at Low Energies 

At low electron energies, the dominant losses will be due 
to Coulomb collisions. Coulomb losses have the property that 
they vary only rather slowly (logarithmically) with 7 (eq. ^]). 
Here, this slow variation is ignored, and the Coulomb losses 
are written as £>(7>0 = bqmii ~ be, where be is a constant given 
by the value in equation (ffl for some appropriate value of 7, say 
7 ~ 10 2 . Then, equation (|l4|) for the change in 7 for a single 
particle gives 

7 = 7;~7w, (32) 
where the characteristic value -ft m , is defined by 



7/o.r = bc&t . 



(33) 



The equation for the evolution of the electron population 
(eq. O) reduces to 



jV(7,0«jV[7 + 7/ 



(34) 



For a short time At <C 7 /be, the population is nearly the initial 
one. However, the Coulomb losses at low electron energies will 
tend to be rapid. Thus, after a time Af 3> 7 /be, the population 
will approach 

jV( 7 ,f)~jVbW,]. (35) 

That is, the population at low energies will vary with time, but 
will be nearly independent of 7, 



N(7) « Ni ow (t) = constant . 



(36) 



3.2. Steady Injection Only 

As an alternative to the case where there is an initial popula- 
tion of particles but no continual source, we consider the case 
where there is continuous injection of particles at a constant rate 
(2(7), but where there is no initial electron population. 

3.2.1. Steady-State Solutions 

At energies where the time scale for losses by particles f/ OSJ 
is much shorter than the age of the cluster, the population will 
tend towards a steady-state distribution. From Figure ||, this is 
likely to occur at high energies (7 > j max ) due to rapid IC and 
synchrotron losses. It may also occur at low energies (7 < ji ow ) 
due to rapid Coulomb losses. If f; 0JJ = [7/^(7)] <C Af, then 
it is likely that dN(j,t)/dt < d[b('j,t)N(j,t)]/dj. In steady- 
state (d/dt = 0), the equation for the evolution of the electron 
population (eq. [5]) becomes 



d", 



[fc(7W(7)] = -G(7)- 



(37) 



If there is no flux of particles from infinite energies and if 
lim 7 ^oo b{-f)N{"f) — * 0, then the steady-state solution is given 
by 

1 f°° 

m ^=TT-J G(7V7'- (38) 

It is useful to consider a source which is a power-law func- 
tion of the electron energy, with 



e(7)=Gi7" 



(39) 



with p > 1 . Then, at high energies where IC and synchrotron 
losses dominate and we can write the loss function as ^(7) « 
&17 2 , the resulting electron distribution is 



7 



-(P+i) . 



(40) 



that is, the resulting electron distribution is one power steeper 
than the source. At low energies where Coulomb losses dom- 
inate, the losses depend only logarithmically on the electron 
energy. If we ignore this slow dependence and take £(7) « be 
(a constant), then the resulting steady-state distribution is 



2i 



bc(p-l) 



7 



<P-r> . 



(41) 



that is, the resulting electron distribution is one power flatter 
than the source. These classical results were derived first by 
Ginzburg (e.g., Ginzburg & Syrovatskii 1964). 

3.2.2. Self-Similar Solutions 

For simple approximations for the loss function and simple 
expressions for the rate of injection Q(j), one can often find 
self-similar solutions which connect the steady-state population 
with non-steady populations which are time dependent. In these 
solutions, the overall normalization of the electron population 
varies with time, but the shape of the energy spectrum remains 
the same if 7 is scaled by some characteristic value. 

At high energies where IC and synchrotron losses dominate, 
one can approximate the loss function by equation ( 16 ), at least 
at low redshifts. The characteristic value of 7 at high ener- 
gies is 7„ rav (eq. [jj]). Then, if the particle injection follows a 
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power-law distribution (eq. 39 ), the self-similar solution for the 
electron population is given by 



N(j,t): 



Qi 



biip-l) 



-(P+D 



1- (1-7/7* 



X) 



1 



7 < In 



(42) 

At low energies where Coulomb losses dominate, the loss 
rate in nearly independent of 7, bhf) w Z?c- Then, the character- 
istic value of the 7 is 7; OM , (eq. ^3[). For a power-law spectrum 
of injected particles (eq. [39]), the self-similar solution for the 
electron population is given by 



#(7,0 = 



bdp-l) 



7 



<P-r> 



l-(l+TW/7) 



-(p-1) 



(43) 



At energies where losses are not important (7 -C 7,„ nv in 
equation |£|] and 7 ^> ji ow in equation [Q), the two self- 
similar solutions are identical, N(<y,t) w Q{"f)t. That is, the 
population is just the sum of the injected particles at that en- 
ergy, as expected if losses have not affected the energies of the 
particles. If ji ow <C 7 mflj: , then the two solutions can simply be 
joined in this central region. With some algebra, one finds that 
this solution can be written as 



#(7,0 ■ 



Qa~ p t ( n 



1 



7 



7 

^fmax 



7/o 



-(p-i) ' 



1 + — .(44) 



This solution applies for all 7 as long as 7/ 0H . <C "fmax- 

4. SOURCE FUNCTION FOR ELECTRONS 
4. 1 . Injected Electron Spectrum 

Models will be calculated which include an initial popula- 
tions of relativistic electrons and/or continuous injection of rel- 
ativistic electrons during the history of the cluster. Although the 
models don't depend on the details of the processes which ac- 
celerate relativistic electrons in clusters, the basic assumption 
is that these particles result from shock or turbulent accelera- 
tion. The shocks and turbulence might be associated with the 
hydrodynamical evolution of the ICM; for example, the initial 
population might be due to an accretion shock during the forma- 
tion of the cluster, and semi-continuous particle injection might 
result from subcluster mergers. Alternatively, the shocks might 
be due to the expansion of radio sources within the ICM. 

The condition for the shocks due to the hydrodynamics of the 
intracluster medium are similar to condition in Galactic super- 
nova remnants (SNRs), except that the length scales are much 
larger and the densities are smaller. Over a significant range of 
relativistic energies, there is considerable theoretical and obser- 
vational evidence that such shock acceleration produces parti- 
cles with a power-law distribution (e.g., Blandford & Eichler 
1987; Jones & Ellison 1991). Thus, we will consider models in 
which the initial electron spectrum at t = ti is a single power-law, 



For typical strong shock compressions, one expects values of 
p = pqk, 2.3 (e.g., Blandford & Eichler 1987; Jones & Ellison 
1991). 

One complication in applying the standard results on shock 
acceleration to the ICM is that we need to determine the popu- 
lation of relatively low energy cosmic ray electrons (7 ~ 300) in 
order to calculate their EUV emission. Such particles have ki- 
netic energies (^150 MeV) for which electrons are relativistic, 
but protons and other ions are not. There are a number of simple 
arguments which suggest that the energy spectrum of superther- 
mal but subrelativistic particles is a flatter power-law than for 
fully relativistic particles (Bell 1978b). The suggestion is that 
the momentum spectrum of particles accelerated in shocks is the 
same power-law for both relativistic and nonrelativistic parti- 
cles and for both ions and electrons. That is, Q(P)dP oc P~ p dP, 
where P is the particle momentum and Q(P)dP gives number 
of particles accelerated with momenta between P and P+dP 
per unit time. For highly relativistic particles E = Pc, but for 
non-relativistic particles the kinetic energy E = P 2 /2m, where 
m is their mass. Under the same assumptions which result in 
a shock acceleration energy spectrum of highly relativistic par- 
ticles with Q(E) oc E~'\ the energy spectrum for superthermal 
but subrelativistic particles is Q(E) oc E~ < - p+1 ^ 2 . For a typical 
relativistic injection spectrum with p w 2.3, the subrelativistic 
spectrum would have a lower exponent pi w 1.65. At energies 
E ~ 100 MeV, this would only affect the shock acceleration 
spectrum of ions. Bell (1978b) suggests that would result in 
the total number of accelerated electrons and ions (at energies 
above m e c 2 ) being nearly equal, while the number of ions with 
E 3> m p c 2 would be larger than that of electrons at the same 
energies. 

However, it is also possible the the spectrum of shock ac- 
celerated electrons would also be affected. Thus, we will also 
consider models in which the initial electron spectrum at t = t\ 
is a broken power-law, 



N(-y,ti) = Na~ 



(45) 



and/or the continuous source of injected electrons is a power- 
law, 

e(7,0 = 2i(f)7~"- (46) 



7 



7 < Ibr 



N( 1 ,t i )=N l x{ (47) 

7>7* r 

and/or the continuous source of injected electrons is a broken 
power-law, with 



7 



2(7,0 = Gi(0 



lb, 



(if 



7 < Jbr 



7 > Jbr 



(48) 



4.2. Normalization of the Electron Population 

Equation || for the evolution of the population is linear in 
#(7), so the solution can be rescaled by increasing or decreas- 
ing N(j,t), the initial value N(j,ti), and the electron injection 
rate (2(7, by any constant factor. Thus, the normalization of 
the solution is arbitrary. I have chosen an arbitrary normaliza- 
tion; however, it is at least somewhat consistent with energetic 
and stability arguments concerning the ICM medium, with the 
observed EUV fluxes of clusters, and with the observed prop- 
erties of shock acceleration in SNRs in our Galaxy. If the ob- 
served EUV radiation from clusters results from IC scattering 
of CMB photons by cosmic ray electrons, the observed EUV lu- 
minosities imply that clusters contain cosmic ray electrons with 
7 J> 300 which have a total energy of the order of £cR,e ~ 10 62 
ergs. It is also unlikely that the cosmic ray energy significantly 
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exceeds the thermal energy of the ICM, since this would likely 
result in the ICM being expelled from the clusters. Rich clusters 
contain approximately 1O 14 M0 of hot gas at a typical tempera- 
ture of 7 x 10 7 K, which implies that the thermal energy content 
of the gas is approximately E gas ~4 x 10 63 ergs. Thus, if the 
cosmic ray electrons had a total kinetic energy of Ecu, e ~ 10 62 
ergs, this would represent approximately 3% of the thermal en- 
ergy content of the ICM. Even when one includes the contribu- 
tion from cosmic ray ions, this ratio is likely to remain of order 
unity or smaller. Finally, we note that the radio observations of 
Galactic SNR imply that shock acceleration produces cosmic 
ray electrons which contain at least a few per cent of the shock 
energy (e.g., Blandford & Eichler 1987). Since essentially all 
of the thermal energy content of the ICM is the result of shock 
heating (§ [j] above), these shocks should also have accelerated 
relativistic electrons containing at least a few per cent of the 
thermal energy content of the intracluster medium. 

I have chosen to normalize the electron spectra so that the to- 
tal amount of kinetic energy injected in electrons with 7 > 300 
is £^ f = 10 63 ergs. In the models which seem most likely to 
represent real, present day clusters, this leads to present day 
electron populations with EcR, e ~ 10 62 ergs. 

4.3. Self-Similar Cluster Accretion Shock 

In most of the models with continual injection of particles, 
we assume for simplicity that the injection rate is constant over 
the time that it acts, dQ(-f,t)/dt = 0. As an alternative, we con- 
sider a simple model for the cosmological evolution of 2(7,0- 
Bertschinger (1985) has given a simple, self-similar solution 
for the secondary accretion of intracluster gas (with or with- 
out associated dark matter) onto a previously collapse cluster 
core. This solution applies only in an Einstein-de Sitter Uni- 
verse, with il = 1 and A = 0. In the Bertschinger solution, the 
accretion shock radius expands with time as r s cx f 8 / 9 , and the 
preshock density varies with the average cosmological density 
p cx f" 2 . The shock velocity obviously varies at v s cx r s /t cx f -1 / 9 . 
In this solution, the preshock gas is cold and the accretion shock 
is always very strong. Thus, rate of shock energy input varies 
as E shock = ( 1 /2) /9 v 3 4 7 rr 2 cx r 5 /9 . 

If we assume that a given fraction of the shock energy goes 
into accelerating the relativistic electrons, then it is useful to 
give the rate of shock energy input at any time in terms of the 
total shock energy over the life time of the cluster. Integration 
of the shock energy equation leads to the following expression 
for the accretion shock energy input in the Bertschinger self- 
similar solution: 

E shock (z) = 4l (1+z) 5 / 6 
E shock 9f l-(l+z,)- 2 / 3 ' 

where t„ is the present age of the Universe, z is the redshift at 
the time when the shock energy is being evaluated, and zi is the 
redshift at which the cluster formed. 

5. NUMERICAL MODELS 

5.1. Techniques 

A number of numerical models for the integrated energy 
spectra of cosmic ray electrons in clusters of galaxies have been 
calculated. The models were constructed so that the present day 
electron populations could be determined up to at least 7 = 10 5 . 
In models without continual electron injection, the present day 
electron populations are generally cut-off at 7„ MV < 10 5 . In the 



models with continual electron injection, the electron popula- 
tion is in steady-state above 7 = 10 5 in all of the models consid- 
ered, and the populations can easily be determined from equa- 
tions ( f38"l ) or (|40|). Also, the rapid loss rates for high energy 
electrons with 7 > 10 s mean that numerical solutions for these 
electrons are time consuming. 

For models with an initial electron population at redshift Zi 
but no subsequent particle injection, the population was cal- 
culated by following the energy losses of individual particles. 
That is, equation (Q) was integrated backward in time to deter- 
mine the initial Lorentz factor 7, corresponding to the present 
day 7 for a large number of particles with energies in the 
range 1 < 7 < 10 5 . The loss equation was integrated using the 
Bulirsch-Stoer technique with adaptive step-size control (e.g., 
Press et al. 1986). The equation of cosmic dynamics (eq. £3|) 
was used to relate redshift and co-moving time. All of the loss 
processes discussed in § were included. Then, the present 
day electron spectrum was determined from particle conserva- 
tion through equation (|l3|). 

For models with continual injection but no initial popula- 
tion, the equation for the evolution of the electron energy spec- 
trum (eq. |j) was integrated with a finite-difference technique 
on a fixed grid of energies covering the range 1 < 7 < 10 5 . I 
found that the energy spectrum at very high energies was well- 
represented by the self-similar solution (eq. ^2|). I assumed this 
solution applied at high energies beyond the upper end of the 
grid, and it was used to determine the flux of particles into the 
grid from higher energies. The time step was regulated so that 
the flux of particles into each grid cell from higher energies was 
a small fraction (< 1 %) of the population of that cell, and so that 
the redshift step was also very small. 

Because the equation for the evolution of the electron popula- 
tion is linear (eq. gj), models with both an initial populations and 
subsequent particle injection were created by the superposition 
of solutions with only and initial population and only subse- 
quent injection. 

Unless otherwise noted, the average thermal electron density 
in the ICM was taken to be n e = 0.001 cm" 3 , the Hubble con- 
stant was Ho = 65 km s" 1 Mpc -1 , and the magnetic field was 
B = 0.3 /iG. In all of the models, the deceleration parameter 
was qo = 0.5. 

In all of the models, there is a population of particles with 
7 w 1, which represent initially relativistic electrons which have 
returned to the thermal population through energy losses. Since 
the number of cosmic ray electrons is small compared to the 
number of thermal electrons, this population is probably not 
important. 

Some of the assumed properties of the numerical models are 
summarized in Table [l]. The models start at a redshift of Zi- In 
models with continual particle injection, this starts at z s - The 
value of the Hubble constant in the models is given in column 
4. The average thermal electron density n e and magnetic field 
B are given in columns 5 and 6. A "Y" in the columns labeled 
"Initial" or "Cont." indicates whether the model includes an ini- 
tial population of relativistic electron at Zi and/or whether there 
is continual particle injection since z s . The values of po, p, j^r, 
and pi indicate the form of the power-law or broken power- 
law injection spectrum, according to equations (p3|)-(p^. If no 
value is given for 7;, f , the injection spectrum is a single power- 
law (eqs. p3]-jl6|). For models with both an initial population 
and continual particle injection, the value of Fj„j is the fraction 
of the total particle energy which is contributed by the contin- 
ual particle injection. The last column gives a comment on the 
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Table 1 

Numerical Models for the Electron Spectrum 



Model 


Zi 


z s 


Ho 

(km s -1 MpcT 1 ) 


n e 
(cm 3 ) 


B 

(pG) 


Initial 
Pop. 


Po 


Cont. 
Inj. 


P 




Pi 


F'm j 

(%) 


Comment 


1 


2.00 


2.00 


65 


0.001 


0.3 


N 




Y 


2.3 






100 


Steady injection 


2 


1.00 


1.00 


65 


0.001 


0.3 


N 




Y 


2.3 






100 


Steady injection 


3 


0.50 


0.50 


65 


0.001 


0.3 


N 




Y 


2.3 






100 


Steady injection 


4 


0.30 


0.30 


65 


0.001 


0.3 


N 




Y 


2.3 






100 


Steady injection 


5 


0.10 


0.10 


65 


0.001 


0.3 


N 




Y 


2.3 






100 


Steady injection 


6 


0.01 


0.01 


65 


0.001 


0.3 


N 




Y 


2.3 






100 


Steady injection 


7 


2.00 


2.00 


65 


0.0001 


0.3 


N 




Y 


2.3 






100 


Model 1 , but lower density 


8 


2.00 


2.00 


65 


0.000 


5.0 


N 




Y 


2.3 






100 


Model 1 , but stronger magnetic field 


9 


0.01 




65 


0.001 


0.3 


Y 


2.3 


N 











Initial pop. only 


10 


0.10 




65 


0.001 


0.3 


Y 


2.3 


N 











Initial pop. only 


1 1 


0.30 




65 


0.001 


0.3 


Y 


2.3 


N 











Initial pop. only 


12 


0.50 




65 


0.001 


0.3 


Y 


2.3 


N 











Initial pop. only 


13 


1.00 




65 


0.001 


0.3 


Y 


2.3 


N 











Initial pop. only 


14 


0.30 




65 


0.0001 


0.3 


Y 


2.3 


N 











Model 11, but lower density 


15 


0.30 




100 


0.001 


0.3 


Y 


2.3 


N 











Model 1 1, but larger H Q 


16 


0.30 




65 


0.001 


0.0 


Y 


2.3 


N 











Model 11, but no magnetic field 


17 


0.30 




65 


0.001 


1.0 


Y 


2.3 


N 











Model 11, but stronger magnetic field 


1 5 


n in 
U. 5U 




CO 




J.U 


V 

1 


Z.J 


M 

IN 








u 


Model 1 1 , but stronger magnetic field 


19 


0.30 




65 


0.001 


5.0 


Y 


2.3 


N 











Model 11, but stronger magnetic field 


20 


2.00 


2.00 


65 


0.001 


0.3 


N 




Y 


2.3 






100 


Bertschinger model 


21 


2.00 


2.00 


65 


0.001 


0.3 


N 




Y 


2.3 


2000 


1.3 


100 


Broken power-law injection 


22 


0.30 




65 


0.001 


0.3 


Y 


2.3 


N 




2000 


1.3 





Broken power-law initial pop. 


23 


0.05 


0.05 


65 


0.001 


0.3 


N 




Y 


2.3 






100 


Steady injection 


24 


0.30 


0.05 


65 


0.001 


0.3 


Y 


2.3 


Y 


2.3 






50 


Both init. pop. & current inj. 


25 


0.30 


0.05 


65 


0.001 


0.3 


Y 


2.3 


Y 


2.3 






25 


Both init. pop. & current inj. 


26 


0.30 


0.05 


65 


0.001 


0.3 


Y 


2.3 


Y 


2.3 






10 


Both init. pop. & current inj. 


27 


0.30 


0.05 


65 


0.001 


0.3 


Y 


2.3 


Y 


2.3 






1 


Both init. pop. & current inj. 


28 


0.30 


0.05 


65 


0.001 


0.3 


Y 


2.3 


Y 


2.3 






0.1 


Both init. pop. & current inj. 



model. 



5.2. Solutions 



Steady Injection 
No Initial Pop. 
Power Law p = 2.3 
z, = 3.00 




FIG. 4. — The relativistic electron population in a cluster in which particles 
have been injected at a constant rate since zt = 2. The injected particles had a 
power-law distribution with p = 2.3. The solid curve (Model 1) is the distri- 
bution in a cluster with an average thermal electron density of n e = 10~ 3 cm -3 . 
The two short-dashed lines are power-law distributions with p' = 1.3 and 3.3, 
which match the resulting particle distributions at low and high energy, respec- 
tively. The two long-dashed curves are similar to Model 1 , but have a decreased 
electron density of n e = 1CT 4 cm -3 (Model 7) or an increased magnetic field of 
B = 5 p,G (Model 8). 



Some of the results for the numerical models for the electron 
population in clusters are given in Table [| For each model, 
the second column gives the total number of relativistic elec- 
trons N, , at present (at z = 0). Columns 3 and 4 give the kinetic 
energy in relativistic electrons at present (z = 0), either for all 
of the electrons (7 > 1) or only for those with 7 > 300. The 
fifth column gives the energy-averaged mean value of 7, de- 
fined here as 

_ J>(7)7 2 a?7 



(7) 



/ ^(7)7^7 



(50) 



Sarazin & Lieu (1998) show that this quantity is useful for es- 
timating the emission properties of these electro ns. Fo r models 
without particle injection at the present time (§ 5.2.2), the last 
column gives the upper cut off to the electron distribution, 7 mar . 

5.2. 1 . Steady Particle Injection 

I first considered models in which relativistic electrons are 
injected into the ICM at a constant rate over the history of the 
cluster, and there was no initial population of electrons. Fig- 
ure Q shows the result of steady injection since an initial red- 
shift of Zi = 2; the solid curve is for Model 1 in Table |lj For 
this large redshift, the age of the cluster exceeds the loss time 
of the electrons for all values of 7 (Figure ||). This is particu- 
larly true when one includes the increase in the energy density 
of the CMB with redshift. For example, Figure [3] shows that IC 
losses will remove all of the particles produced at z ~ 1 down to 
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low values of 7 where Coulomb losses dominate and the life- 
time is short. Thus, one expects that the particle distribution 
will approach steady state, with the rate of injection balancing 
the rate of loss at most energies. The injected electrons had a 
simple power-law distribution (equation^) with p = 2.3. Fig- 
ure l] shows that the resulting electron population is essentially 
a broken power-law. The dashed lines show power-law dis- 
tributions with indices of p' = 1.3 and 3.3. At lower energies 
where Coulomb losses dominate (7 < 100), the distribution is 
very nearly a power-law with p' =1.3. This is one power flatter 
than the injected population, as expected from the steady-state 
result of equation (ffTl). At high energies where IC and syn- 
chrotron losses dominate (7 ;> 1000), the distribution is nearly 
a power-law with p' = 3.3. This is one power steeper than the 
injected population, as expected from the steady-state result of 
equation (Q. 



Table 2 

Results for the Electron Energy Spectrum 



EcR.e 

7 > 1 7 > 300 



7». 



Model 


(10 66 ) 


(10 6 


ergs) 






1 


11.98 


28.39 


9.48 


544.6 




2 


14.96 


35.46 


11.85 


544.6 




3 


21.16 


50.15 


16.81 


546.2 




4 


29.12 


67.32 


23.50 


564.2 




5 


65.38 


119.23 


42.59 


675.7 




6 


377.34 


264.52 


71.55 


1159.4 




7 


95.67 


91.58 


11.50 


204.7 




8 


11.73 


21.10 


4.29 


302.0 




9 


118.87 


198.19 


62.98 


761.9 


15107 


10 


7.67 


56.76 


23.17 


334.0 


1379 


11 


1.95 


13.38 


0.70 


153.6 


338 


12 


0.67 


2.17 


0.00 


60.6 


110 


13 


0.00 


0.00 


0.00 


3.1 


83 


14 


44.38 


88.65 


4.17 


93.2 


405 


15 


3.93 


28.99 


7.67 


209.9 


582 


16 


1.95 


13.43 


0.74 


154.4 


340 


17 


1.93 


12.82 


0.31 


147.0 


320 


18 


1.76 


9.01 


0.00 


104.7 


211 


19 


1.44 


4.66 


0.00 


61.6 


114 


20 


8.51 


21.01 


6.93 


527.5 




21 


0.54 


11.12 


8.13 


1513.9 




22 


0.61 


7.57 


0.93 


212.0 


333 


23 


113.10 


159.80 


53.50 


787.0 




24 


57.52 


86.59 


27.10 


740.1 




25 


29.74 


49.98 


13.90 


664.4 




26 


13.06 


28.02 


5.98 


521.8 




27 


3.06 


14.84 


1.23 


224.6 




28 


2.06 


13.52 


0.75 


161.5 





The approach to steady-state in models with constant injec- 
tion is shown in Figure which gives the present day pop- 
ulations in models with initial redshifts of z,- = 0.01, 0.1, 0.3, 
0.5, 1, and 2 (top to bottom). (These are Models 1-6 in Ta- 
ble [l].) All of the models have an average electron density of 
n e = 10~ 3 cm" 3 , and a magnetic field of B = 0.3 ^G. The parti- 
cles are injected with a power-law spectrum with p = 2.3. The 
older models (z,- = 0.5, 1, and 2) have achieved steady-state, 
since the loss times of electrons of all energies are less than the 
age of the cluster. The younger cluster models show departures 
from steady-state. For these models, the energy spectra of elec- 
trons can be divided into three regions. At very high energies 



(7 ;> "/ ma x), the IC losses are sufficiently rapid that the loss time 
scale (eqn. [n]) is shorter than the age, and the populations ap- 
proach steady-state, and the population is a power-law which 
is one power steeper than the injected spectrum. Similarly, at 
low energies the loss time scale is also generally shorter than 
the age (Figure fy. The population also approaches steady-state 
below some lower value of 7 <J 7; 0M „ and population is a power- 
law which is one power flatter than the injection spectrum. At 
intermediate energies ji mv <C 7 <C "fma.i, the loss time scale is 
longer than the age, and the population is given by the accu- 
mulation of the injection rate. Thus, it has the same spectrum 
as the injection spectrum. For younger injection models where 
"flow <C Jmax> the resulting electron population isgiven quite ac- 
curately by the self-similar solution (equation p4|). 




log 7 



FIG. 5. — The present day relativistic electron populations in a series of 
models with steady particle injection but no initial population of particles. The 
solid curves show models for clusters which started at redshifts of Ziff 2, 1, 0.5, 
0.3, 0.1, and 0.01 (bottom to top). These are Models 1-6 in Table M. In all of 
the models, the particles are injected at a constant rate with a power"aw distri- 
bution with p = 2.3. The short-dashed curve gives the total power-law spectrum 
of all of the injected particle over the cluster lifetime. 



Figure [| also shows the effect of changing the gas density 
or magnetic field in the models. The two long-dash curves are 
the same as Model 1, but with the electron density decreased 
to n e = 10~ 4 cm" 3 (Model 7) or the magnetic field increased to 
B = 5 /iG (Model 8). Increasing or decreasing the gas density 
mainly increases or decreases the rate of Coulomb losses at low 
energies, and this decreases or increases the electron popula- 
tion at low energies. Reducing n e also decreases the electron 
energy or value of 7 at which the transition between the two 
power-law slopes occurs in steady-state. Similarly, increasing 
magnetic field enhances the rate of synchrotron losses at high 
energies. This reduces the population of at high energies, but 
the effect is only important if the field exceeds about 3 /iG. 

5.2.2. Initial Population with No Later Particle Injection 

Next, models were calculated with an initial electron popu- 
lation at a redshift z,-, but with no subsequent injection of parti- 
cles. Figure ^| shows the result for models with differing values 
of the initial redshift z,-. The initial particle population was a 
single power-law with po = 2.3 (eqn. [45]) At high energies, the 
electron population is reduced and the shape of the spectrum 
steepened by IC and synchrotron losses. There is an upper cut- 
off to the electron distribution (eqns. |l9|-|20|) for values of 7 
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beyond a maximum values j max , given approximately by equa- 
tion (27 ), or equation ( 18 ) for low redshifts. The cutoff energies 
in the numerical models are in reasonable agreement with these 
analytic approximations as long as ^ max ^> 300. The value of 
7„ lflA decreases as the age of the electron populations (or, equiv- 
alently, the initial redshift zi) increases. 



Initial Pop. 
No Steady Injection 
Power Law p = 2.3 
n e = 10 cm 
= 0.3 fjG 




log 7 



FIG. 6. — The relativistic electron population in models with an initial pop- 
ulation of particles generated at a redshift Zh Dut no subsequent injection. The 
initial population, which is shown as a dashed line, had a power-law distribu- 
tion with po = 2.3. The solid curve show the resulting present population for 
values of the initial redshift of z, = 0.01, 0.1, 0.3, and 0.5. 



spectrum at low energies (eqn. [3(|). All of the older cluster 
models in Figure || show this behavior. In the younger models 
(Zi <, 0.3), there is an extended region 7/ ovv <J 7 < j max where the 
initial population is preserved. In the older models (z, J> 0.3), 
the flat energy spectrum at low energies and steep cut-off at 
high energies nearly meet, and the energy spectra have a "top 
hat" form. 

Figure ^ shows the effect of varying two of the basic pa- 
rameters of the models. The solid curve is Model 11, with an 
average electron density of n e = 10~ 3 cm" 3 and a Hubble con- 
stant of Hq = 65 km s -1 Mpc -1 . In all of these models, the 
electron population was injected into the cluster at z,- = 0.3. The 
short-dashed curve (Model 14) is a model with all of the same 
parameters, but with a lower value for the electron density of 
n e = 10~ 4 cm" 3 . Reducing the electron density lowers the rate 
of Coulomb losses, which dominate at low energies. As a re- 
sult, increasing or decreasing the electron density decreases or 
increases the population of low energy electrons. 

Changing the value of the Hubble constant changes the age 
of the cluster for a fixed value of z, and qo. The long-dashed 
curve in Figure [7] shows Model 15, which has all of the same 
parameter, except the Hubble constant is increased to Ho =100 
km s -1 Mpc" 1 . Increasing the Hubble constant shortens the age, 
and reduces the effect of losses. Thus, the electron population 
is larger at all energies, and the cutoff value of 7 mflJC is increased. 
Lowering the Hubble constant has the opposite effect. 



In order to have any significant population of primary elec- 
trons with 7 > 10 2 at the present time, Figures and @ show 
that there must have been a substantial injection of particles 
into clusters at moderately low redshifts, z <J 1 . 



n. — 10 cm 



Initial Pop. 

No Steady Injection 

Power Law p„ = 2.3 

z, = 0.3 

B = 0.3 /zG 




FIG. 7. — The effect on the electron population of changing several of the 
parameters. The solid curve is Model 1 1 for an initial population of particles 
generated at a redshift zi = 0.3, but no subsequent injection. The short-dashed 
curve (Model 14) shows the result of reducing the average electron density 
from 10~ 3 to 10 -4 cm . The long-dashed curve (Model 15) shows the result of 
increasing the Hubble constant from 65 to 100 km s _1 Mpc -1 . 



At low energies (7 <^ 7/ OH ), the electron population is reduced 
and the shape of the distribution flattened by C oulomb losses to 
the thermal plasma. As discussed in § 3.1.6, the slow varia- 



tion of the Coulomb loss function with 7 results in a flat energy 




Initial Pop. 

No Steady Injection 

Power Law p = 2.3 

z, = 0.3 

n e = 10 cm 



1.5 
log y 



FIG. 8. — Models for relativistic electrons with a variety of intracluster mag- 
netic field strengths (Models 11, 16-19). The curves show the populations for 
magnetic fields of 0, 0.3, 1, 3, and 5 fiG. The curves for no magnetic field and 
0.3 fj,G are nearly indistinguishable. 



I also constructed models with a variety of average values 
for the intracluster magnetic field, which can affect the electron 
population through synchrotron losses. Figure [| shows Mod- 
els 16, 11, and 17-19, with magnetic fields of B = 0, 0.3, 1.0 
3.0, and 5.0 ^G. All of the other parameters of Models 16-19 
are identical to Model 1 1 . Because synchrotron losses have the 
same variation with 7 as IC losses, they affect the electron pop- 
ulation in a very similar way. Unless the ICM magnetic field 
exceeds 1 /iG, synchrotron losses do not have a very signifi- 
cant effect. For larger values of B, the number of high energy 
electrons and the value of j max are reduced as B is increased. 
The values of 7,„ m - are in good agreement with those predicted 
by equation (bTJ). However, to produce an appreciable effect, 
the magnetic field is required to be considerably larger than 
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expected in clusters (Kim et al. 1990, 1991: Kronberg 1994; 



Fusco-Femiano et al. 1998). 
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FIG. 9. — The solid curve is the relativistic electron population in a model 
based on the iiertschinger (1985) self-similar solution for secondary infall in a 
cluster (eq. [tt9p, Model 20). For comparison, the dashed curve is the steady 
injection model with the same properties (Model 1). 



5.2.3. Self-Similar Cluster Accretion Shock 

I also calculated a model in which particles were injected 
with a rate proportional by the rate of deposition of shock 
energy (eq. ^9|) in the self-similar secondary infall model of 
Bertschinger (1985). In Figure ^, this model (Model 20) is 
compared to the steady injection model (Model 1) in which 
the rate of injection of particles is a constant. The two mod- 
els give very similar results. The small differences between 
the self-similar secondary infall model and the steady injec- 
tio n m odel are mainly due to differences in the normalization 
(§ 4.2). The models are normalized based on the total input 
of relativistic electrons, and the Bertschinger model and the 
steady model have different variations in the rate of injection 
with time or redshift (eq. ^). If the present day rate of injection 
is the same, the maximum difference between the Bertschinger 
model and the steady injection model is only about 10%, and 
occurs near the knee in the particle spectrum (at 7 « 160 in 
Figure ^). The very near agreement of the Bertschinger model 
and the steady injection model reflects the fact that both mod- 
els have nearly reached steady-state, and the current electron 
spectrum is mainly determined by the present rate of injection, 
rather that the past rate over cosmological time scales. 



5.2.4. Models with Broken Power Law Particle Injection 



In § 4.1 , we noted that the spectrum of cosmic ray protons 
produced by shock acceleration is expected to flatten some- 
what at kinetic energies E <; 1 GeV, where the protons become 
nonrelativistic. In most models, it is expected that the elec- 
tron spectrum will remain nearly a power-law to much lower 
energies. Here, we consider that possibility that the spectrum 
of electrons also flattens below 1 GeV. I adapt the simple two 
power-law form of equations [|7]) and (f48|). In Figure [l(J the re- 
sulting electron populations in models with broken power-law 
models for the initial population or for the injection spectrum. 
For models with steady injection (Model 21), the broken spec- 
trum produces a low energy electron spectrum which is a flatter 
power-law by the same amount that the injection spectrum if 



flattened at low energies. For models with initial broken power- 
law spectra, the population is very flat at low energies, and has 
the same high energy cut-off as the model with a single power- 
law. However, the broken power-law models have a peak at 
energies just below j max if ^ br > % nax . 




log 7 



FIG. UL — Models with broken power-law injection (eq. ftcj) or initial spec- 
tra (eq. W7p are compared to models with simple power-laws. The lower curves 
compare a model with an initial population (at n = 0.3 ) which is a broken 
power-law (solid curve, Model 22) with an otherwise identical model with a 
single power-law (dashed curve, Model 11). In these models, there is no sub- 
sequent particle injection. The upper curves compare models with no initial 
population but a steady rate of particle injection since z, = 2. The solid curve 
has a broken power-law injection spectrum (Model 21), while the dashed curve 
is for a single power-law (Model 1). For clarity of presentation, the steady 
injection curves were displaced upwards by one order of magnitude. 

5.2.5. Models with Both an Initial Population and Steady 

Injection 

Finally, we consider models in which there is an initial pop- 
ulation in the cluster at z,-, and in which there also is current in- 
jection of new particles. I assume that this current injection has 
occurred at a constant rate since redshift z s - These two com- 
ponents might represent particle populations which were pro- 
duced by shocks associated with the formation of the cluster 
at Zi and by a current merger event, which started at z s . Be- 
cause the equation for the evolution of the electron population 
is linear (eq. ||), these models are just a linear combination of 
the models discussed in §§ 5.2.1 and |5.2.2| . An illustrative set 
of such models are shown in Figure |1 1|. I assumed z, = 0.3, 
Zj = 0.05, n e = 10~ 3 cm" 3 , and B = 0.3 /iG in these models. Both 



the initial population and current injection are given by a single 
power-law spectrum with p = po = 2.3. The different models are 
characterized by differing values of F/„j, which is the fraction 
of the particle energy which has been injected by the current 
steady injection process, as opposed to the initial population. 
Models are shown with values of F inj = 100%, 50%, 25%, 10%, 
1%, 0.1%, and 0%. (Models 23-28 and 1 1). 

For models with intermediate values of the energy fraction 
due to the ongoing injection, the energy spectrum at low ener- 
gies is flatter than the p' = 1.3 power-law due to steady state 
injection, but steeper then the flat distribution expected for the 
initial population. There is a rapid drop-off in the electron pop- 
ulation at -f max . Above 7 mflx , the electrons are only due to cur- 
rent injection. Of course, in a real cluster it is unlikely that the 
initial population would have been produced at a single redshift. 
Thus, the fall off at -f max is likely to be more gradual, and there 
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might be a more continuous transition to the current injection 
population above j max . However, these models should give a 
general sense of the electron energy spectra to be expected in 
real clusters. 




Initial Pop. (z, = 0-3) & 
Steady Injection (z s — 0.05 
Power Law p = 2.3 
n e = 10 cm 
0.3 liG 



log y 



FIG. 1 1 . — The electron energy spectra in models with both an initial elec- 
tron population starting at zt = 0.3 and continual injection of new particles since 
Z s = 0.05. All of the models have a single power-law spectrum for particle in- 
jection with p = pq = 2.3, an electron density of n e = 10~ 3 cirr 3 , and a magnetic 
field of B = 0.3 pG. The values of the fraction of particle energy due to current 
injection are (top to bottom) F in j = 100%, 50%, 25%, 10%, 1%, 0.1%, and 0%. 



6. INVERSE COMPTON EMISSION 

I also calculated the emission produced by the relativistic 
electron populations due to the inverse Compton scattering of 
CMB photons. Let L v dv be the luminosity of IC emission at 
frequencies from v to v + dv. Then, the spectrum of IC emis- 
sion is related to the electron spectrum N(j) by 



oo pi 

N(j)d~/ / J 
Jo 



where 



J-(x) = \+x+2x\nx-2x 



T(x)dx, (51) 



(52) 



Here, J(v) is the mean intensity at a frequency v of the radiation 
field being scattered. For the CMB, this is just the black body 
function 



J{v)=B v {T CMB )-. 



1 



2hv i 



c 1 exp(hv/kT C MB)-l 



(53) 



Figure [T2| shows the IC spectra of a series of models with a 
steady rate of particle injection, but no initial electron popula- 
tion. These are the same models (Models 1-6) whose electron 
populations were shown in Figure ||. The models shown have 
had steady injection starting at redshifts z; = 2, 1, 0.5, 0.3, 0.1, 
and 0.01 (bottom to top) and continuing to the present. The 
dashed line is the spectrum produced by a model in which the 
electron population is just that injected into the cluster; since 
this involves no evolution, this model is equivalent to z, = 0. 
For this model with no evolution, the spectrum is a power-law 
L v oc v a , where a = — (p— l)/2. This gives a = -0.65 for the 
injection spectrum of p = 2.3. At frequencies above 10 17 Hz, 
the spectrum quickly steepens with time due to IC losses, and 
becomes a half power steeper, a « — (p /2) -1 . 15. At low fre- 
quencies, the spectrum flattens more gradually due to Coulomb 



losses, approaching a spectrum which is a half power flatter, 
a « -{p-2)/2 m -0.15. If the steady injection has occurred 
since a redshift of Zi J> 0.2, the spectrum is near that of a steady- 
state electron population, with these two power-laws meeting at 
a knee at z/- 3 x 10 16 Hz. 




FIG. 12. — The IC emission spectra for cluster models with steady particle 
injectioji but no initial population of particles. The notation is similar to that in 
Figure n. The solid curves show spectra for clusters which started at redshifts 
of zi =p5. 1, 0.5, 0.3, 0.1, and 0.01 (bottom to top). These are Models 1-6 in 
Table Ml. The dashed curve gives the spectrum for z, =0 (i.e., no evolution of 
the electron population due to losses). 




FIG. 13. — The IC emission spectra for cluster models with an initial popu- 
lation of electrons but no subsequent injection. The notation is similar to that in 
Figure H. The solid curves show spectra for clusters which started at redshifts 
of Zi = 0.01, 0.1, 0.3, and 0.5. These are Models 9-12 in Table The dashed 
curve gives the spectrum for z.i = (i.e., no evolution of the electron population 
due to losses). 



The IC spectra of models with an initial electron population 
but no subsequent acceleration are shown in Figure O. This 
Figure gives the spectra for the same models shown in Figure 0, 
in which the initial population of electrons occurs at redshifts 
of n = 0, 0.01, 0.1, 0.3, and 0.5. The electron energy distri- 
butions of these models have a cut off at high energies, 7 mflx , 
which results from rapid IC and synchrotron losses by high en- 
ergy electrons (see eqs. 18, 27], [3l], and Figure ||). As a result, 
the IC spectra have a very rapid fall off at high energies. If 
the initial electron spectrum is a power-law with an index of po 
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(eq. ^5(), and if 7„, cu is large enough that the energy losses are 
dominated by IC and synchrotron losses, then the electron en- 
ergy spectrum at a later time is given by equation (]20|). For this 
electron spectrum, the IC radiation spectrum at high frequen- 
cies, v > 4jl ax kT CM B/h, is given by: 



24irha T ( kT CM B 

po+2 



N\ (2j lnax ) 



-po+l 



hv 



4 lLx kT CMB 



exp 



hv 



4lLx kT CMB 



(54) 



where V is the gamma function. Thus, the spectrum falls off 
exponentially or slightly faster at high frequencies. 

At low frequencies, the IC spectra of these initial popula- 
tion only models flatten and become slowly increasing with 
frequency as the models age. We expect that the low energy 
electron spectrum in an older model with only an initial elec- 
tron population will become nearly independent of 7 (eq. 153 
and Figure^), except for a slowly varying logarithmic factor. If 
N(j) « Ni ow at low energies, where Ni ow is a constant, then the 
IC spectrum is given by 



^^C(5/2)f^ 
5c l \ h 



5/2 



1/2 



(55) 



where C(5/2) « 1 .3415 is the Riemann zeta function. Thus, the 
low frequency IC spectrum varies as L v oc u 1 ' 2 for older clusters 
with initial electron populations only. 
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FIG. 14. — The IC emission spectra for cluster models with both an ini- 
tial population of electrons (at z, = 0.3) and subsequent steady injection (since 
Zs = 0.05). The models are the same as those shown in Figure llll The values 
of the fraction of the particle energy which is due to the current injection are 
(top to bottom) Fi„j = 100%, 50%, 25%, 10%, 1%, 0.1%, and 0%. These are 
Models 23-28 and 11 in Table [l]. 

The spectra of models with both an initial population of par- 
ticles and particle injection at the present time are shown in Fig- 
ure |l4| These are the models whose electron energy spectra are 
given in Figure [fl] (Models 23-28 and 1 1). The different curves 
represent models having differing fractions Fj n j of the total elec- 
tron energy due to the present rate of particle injection. The 
spectra in Figure [l4] correspond to values of Fi n j = 100%, 50%, 
25%, 10%, 1%, 0.1%, and 0%. Obviously, the spectrum of the 
models with either a pure initial electron population (F/„j = 0%) 



or no initial electron population (Fi„j = 100%) are like those 
given in Figures [jj] and [j~2[ respectively. Models in which the 
current rate of particle injection provides a small but signifi- 
cant fraction of the total electron energy have hybrid spectra. 
At low frequencies (v <; 10 17 Hz), they have an extended hump 
of emission, with a rapid fall off above v ~ 10 15 Hz. How- 
ever, they also have an extended hard tail of emission at high 
frequencies, which has a power-law spectrum with a spectral 
index ofa« —(p/2) « -1.15. 
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FIG. 15. — The IC emission spectra for different types of cluster models. 
Models 1 has continual particle injection but no initial electron population (Ta- 
ble hj). Model 11 has an initial electron population, but no subsequent electron 
injection. Model 22 is identical to Model 1 1 , but the initial electron energy 
spectrum has a break at an energy of 1 GeV. Model 27 has both an initial elec- 
tron population, and a source of electrons at present. The scale near the bottom 
of the Figure shows the portion of the electromagnetic spectrum. 



The results for different types of models for the cluster elec- 
tron population are summarized in Figure [l5|. Models 1 has 
continual particle injection in steady state with losses. The ini- 
tial electron population in Model 1 1 was created at z, = 0.3, but 
there has been no subsequent particle injection. Model 22 is 
identical, but the initial electron energy spectrum had a break 
at an energy of 1 GeV. This causes the IC spectrum to decline 
more rapidly toward low frequencies. Finally, Model 27 has 
both an initial electron population (formed at z, = 0.3) and par- 
ticle injection at present (since z, = 0.05) which contributes 1% 
of the total electron population over the history of the cluster. 
At the bottom of Figure ]l5[ there is a scale which shows the por- 
tion of the electromagnetic spectrum involved. The boundaries 
were taken to be: IR-Opt, 7000 A; Opt-UV, 3300 A; UV-X- 
ray, 0.1 keV; X-ray-7-ray, 100 keV. We discuss several of these 
spectral bands separately below. Because many other processes 
contribute to the gamma-ray emission by these particles, I will 
discuss their observable gamma-ray properties in another paper. 
Some of the properties of the IC emission from the models for 
clusters are listed in Table || The total luminosity Lie of the IC 
emission is given in the second column. The other entries are 
discussed below. 

6.1. EUV and Soft X-ray Emission 

Figure [l^ shows that the fluxes from different models all tend 
to agree at frequencies ~2 x 10 16 Hz, which corresponds to a 
photon energy of ^80 eV. Thus, EUV and very soft X-ray emis- 
sion might be expected to be a universal property of clusters of 
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Inverse Compton and Radio Emission from Models 




Uc 


L v (EUV) 


Leuv 


a EUV 


L v (HXR) 




a HXR 


L„ (radio) 


' ^ I'tfll ii' 


L„ (UV) 


a Opt.UV 


Model 


(10 ergs 1 


(10 27 ergs) 


(10 44 eras s" 1 ) 




(10 25 ergs) 


(10 44 ems s" 1 ) 




(10 33 ergs) 




10 28 ergs) 




1 


18.30 


6.406 


1.9685 


-0.68 


1.82 


2.90 


-1.12 


3.41 


-1.14 


2.10 


-0.28 


2 


22.86 


8.002 


2.4590 


-0.68 


2.27 


3.63 


-1.12 


4.26 


-1.14 


2.63 


-0.28 


3 


32.42 


11.352 


3.4884 


-0.68 


3.23 


5.14 


-1.12 


6.06 


-1.14 


3.70 


-0.28 


4 


44.90 


15.202 


4.7490 


-0.64 


4.52 


7.20 


-1.12 


8.49 


-1.14 


4.78 


-0.30 


5 


93.30 


21.628 


6.9003 


-0.60 


11.04 


17.59 


-1.12 


20.75 


-1.14 


7.91 


-0.39 


6 


268.00 


27.269 


8.5045 


-0.64 


44.94 


74.22 


-0.72 


172.59 


-1.08 


13.63 


-0.59 


7 


23.98 


10.906 


3.0394 


-0.89 


1.84 


2.94 


-1.12 


3.42 


-1.14 


7.30 


-0.54 


8 


8.02 


4.166 


1.1759 


-0.86 


0.58 


0.92 


-1.15 


440.03 


-1.12 


1.90 


-0.33 


9 


209.07 


26.360 


8.2593 


-0.63 


38.87 


62.39 


-0.85 


46.61 


-3.62 


12.31 


-0.53 


10 


26.44 


16.356 


5.2954 


-0.57 


0.00 


0.00 




0.00 




4.25 


-0.14 


1 1 


2.86 


4.788 


1.0884 


-1.47 


0.00 


0.00 




0.00 




1.52 


+0.00 


12 


0.18 


0.004 


0.0005 


-11.93 


0.00 


0.00 




0.00 




0.41 


-0.29 


13 


0.00 


0.000 


0.0000 




0.00 


0.00 




0.00 




0.00 




14 


11.59 


14.505 


3.4343 


-1.33 


0.00 


0.00 




0.00 




9.25 


-0.45 


15 


8.48 


10.553 


3.0849 


-0.80 


0.00 


0.00 




0.00 




2.65 


-0.06 


16 


2.89 


4.828 


1.1023 


-1.45 


0.00 


0.00 




0.00 




1.52 


+0.00 


17 


2.63 


4.393 


0.9561 


-1.60 


0.00 


0.00 




0.00 




1.51 


-0.01 


18 


1.31 


1.510 


0.2235 


-3.31 


0.00 


0.00 




0.00 




1.37 


-0.07 


19 


0.40 


0.013 


0.0016 


-11.03 


0.00 


0.00 




0.00 




0.87 


-0.29 


20 


13.18 


4.826 


1.4753 


-0.69 


1.27 


2.03 


-1.12 


2.37 


-1.14 


1.58 


-0.27 


21 


18.32 


2.674 


0.9804 


-0.32 


2.45 


3.90 


-1.12 


4.59 


-1.14 


0.34 


+0.10 


22 


2.23 


4.290 


1.1006 


-1.16 


0.00 


0.00 




0.00 




0.54 


+0.27 


23 


140.48 


24.309 


7.6864 


-0.61 


20.77 


32.97 


-1.10 


39.17 


-1.14 


10.12 


-0.46 


24 


71.65 


14.548 


4.3874 


-0.71 


10.39 


16.48 


-1.10 


19.58 


-1.14 


5.75 


-0.42 


25 


37.26 


9.668 


2.7379 


-0.86 


5.19 


8.24 


-1.10 


9.79 


-1.14 


3.64 


-0.35 


26 


16.62 


6.740 


1.7482 


-1.07 


2.08 


3.30 


-1.10 


3.92 


-1.14 


2.37 


-0.23 


27 


4.23 


4.983 


1.1544 


-1.40 


0.21 


0.33 


-1.10 


0.39 


-1.14 


1.61 


-0.04 


28 


3.00 


4.800 


1.0945 


-1.46 


0.02 


0.030 


-1.10 


0.04 


-1.14 


1.53 


-0.01 



galaxies, if they have all produced populations of relativistic 
electrons with total energies which are at least ~10~ 2 of their 
thermal energies. The one constraint is that at least a portion of 
these particles must have been injected at moderate to low red- 
shifts, Zi <,!■ The fact that essentially all of the models produce 
a strong and similar flux at EUV energies may help to explain 
the fact the excess EUV emission has been detected in all of the 
cluster observed with the Extreme Ultraviolet Explorer (EUVE) 
satellite and which lie in directions of sufficiently low Galactic 
columns that this radiation is observable (Lieu et al. 1996a,b; 
Mittaz, Lieu, & Lockman 1998; Sarazin & Lieu 1998). Val- 
ues of Lv(EUV) at a frequency of v = 2 x 10 15 Hz are listed in 
column 3 of in Table ||. From these values, one finds that the 
emission at EUV energies is fairly directly related to Egg e , the 
total amount of energy injected in electrons with 7 > 300. The 
average relationship is about 

L„(EUV) ~ 6 x 10 27 ergs . (56) 

\ 10 bJ ergs / 

Models which are likely to apply to real clusters have spectrum 
which agree with equation ( pq ) to within a factor of 4. I have 
also calculated the total luminosity in the EUVE observing band 
of photon energies of 65 to 245 eV. These values of Lguv are 
given in column 4 of Table ||. 

While the amount of EUV emission is fairly constant from 
model to model, the spectrum depends strongly on the amount 
of recent particle injection. I have fit the model spectra in the 
EUVE band of 65 to 245 eV to a power-law spectrum of the 
form L v oc v a . Specifically, logL^ was fit to a linear function 
of log v using the least squares method. The best-fit values of 
the energy spectral index aguv are listed in column 5 of Ta- 
ble 0. Large negative values occur for models without current 



particle injection and with small values of 7 ma . t -, where the IC 
spectrum is falling off exponentially (eq. |54|). Thus, a power- 
law may not be a good fit to the spectrum. For most of the 
models, the spectral index lies in the range from about -0.6 to 
-1.5, but very large negative values are also possible. The spec- 
trum drops more rapidly with frequency as the proportion of 
electrons which are currently being injected decreases. This is 
most evident in Models 23-28, which are identical except for a 
decreasing fraction of relativistic electrons due to ongoing par- 
ticle injection. 

Because the EUVE observations of clusters have no spec- 
tral resolution, there is no detailed information on the observed 
spectrum of EUV emission (Lieu et al. 1996a,b; Mittaz, Lieu, 
& Lockman 1998). However, the ratio of EUVE fluxes to 
those in the softest bands of the ROSAT PSPC suggest that the 
EUV spectra are generally steeply declining. The spectra of 
the Coma cluster is less steeply declining than that of A1795; 
this might be understood as the result of particle acceleration by 
merger shocks, since Coma appears to be undergoing a merger 
or mergers (e.g., Burns et al. 1994; Donnelly et al. 1999), while 
A1795 appears regular and relaxed (e.g., Briel & Henry 1996) 
Also, the EUV spectra appear to get steeper with increasing 
radius in the A1795 cluster (Mittaz, Lieu, & Lockman 1998). 
This might be the result of decreasing gas density with radius; 
compare Models 14 and 1 1 in Table |f and Figure 0. 

6.2. Hard X-ray IC Emission and Radio Halos 

At photon energies of 0.5 to 10 keV, the dominant emission in 
most clusters is the thermal emission from the hot ICM, and IC 
emission will be difficult to detect. However, at energies ;>20 
keV, IC emission should again become observable, assuming 
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that it is dropping as power-law function of frequency, while the 
thermal emission drops as an exponential. In Table the hard 
X-ray (HXR) power U at v = 10 19 Hz (41.4 keV, col. 6), the 
luminosity L HX r from 20 to 100 keV, and best-fit HXR power- 
law spectral index chhxr from 20 to 100 keV are given for all of 
the models in Table |l[ IC emission at photon energies of ~50 
keV will be produced by electrons with 7 ~ 10 4 . These par- 
ticles have rather short lifetimes (ti oss 10 9 yr), and are only 
present in clusters in which there has been substantial particle 
injection since z, < 0.1. As Table || shows, only the models 
with current (or very recent, as in Model 9) particle injection 
have any significant HXR emission. 

Because of the short lifetimes of the particles producing 
HXR emission, these electrons are likely to be close to steady- 
state if present in significant quantities. The expected steady- 
state spectral index if IC losses dominated would be ochxr = 
-(p+ 1)/2« — 1.65. The best-fit spectral indices are flatter than 
this, chhxr « — 1.1, mainly because other loss processes are im- 
portant at the lower energy end of the HXR band (Fig. [I]). With 
the exception of the two models in which electrons have only 
been injected in the cluster extremely recently (at z, = 0.01, 
Models 6 and 9), all of the HXR spectral shapes are very simi- 
lar, as expected for steady-state populations. The differences in 
the HXR luminosities just reflect differences in the present rate 
of particle injections. To a good approximation, the present day 
value of Lhxr (20-100 keV) is simply given by 



Lhx«w0.17£ C j?.*(7>5000). 



(57) 



where EcrAj > 5000) is the total present rate of injection of 
energy in cosmic ray electrons with 7 > 5000. The best-fit co- 
efficient (0. 17 in eqn. |5^) depends somewhat on the power-law 
index of the injected electrons; the value of 0.17 applies for 
p = 23. 




log v (Hz) 



FIG. 16. — The radio synchrotron emission from three cluster models. The 
solid curve gives the radio spectrum of Model 23, which is typical of the mod- 
els with current particle injection which has reached steady-state. Almost all of 
the models which have significant radio emission have spectra which are nearly 
identical to Model 23. In Model 6, the electrons haven't reached steady-state, 
and the spectrum is somewhat flatter than that of Model 23. On the other hand, 
in Model 9, all of the electrons were injected at z, = 0.01, and the highest energy 
electrons have lost energy by IC and synchrotron emission. This spectrum is 
much steeper than Model 23 at high frequencies. 



The same relativistic electrons which produce HXR emis- 
sion by IC scattering will also produce radio emission by syn- 
chrotron emission. We have calculated the synchrotron radio 



emission from our models for the intracluster electron popula- 
tion. Averaging over all orientations of the magnetic field and 
all pitch angles for the electrons gives a radio power of 



V3e 3 B , 
Lu = / N(j)R[x(j)]d~/ . 



(58) 



The function R(x) is defined as (Ghisellini, Guilbert, & Svens- 
son 1988) 



R(x) = 2x 2 I K 4/3 (x)K 1/3 (x)- ~x [KfoW-KfoV) 



(59) 



where K is the modified Bessel function. The normalized fre- 
quency variable x is 

(60) 



where v c is the cyclotron frequency 

eB 



2nm e c 



(61) 



The values of the magnetic fields in the models are given in 
Table |l| For an electron with a given energy, the radio spec- 
trum in equation ( |5S| ) is maximum for x « 0.1146. Elec- 
trons with a Lorentz factor of 7 produce radio emission with 
v — 100(B//^G)(7/10 4 ) 2 MHz, and rather high electron ener- 
gies (7 ;> 10 4 ) are needed to produce observable radio emission. 
Thus, these are likely to be the same electrons which produce 
the hardest HXR emission. 

Column 9 of Table | gives the resulting radio powers of the 
models at a wavelength of 91 cm (y = 328 MHz), while col- 
umn 10 gives the spectral index between 91 and 21 cm (v = 1 .4 
GHz). The results are very similar to those for the HXR emis- 
sion. Because of the short lifetimes of the high energy electrons 
required for radio emission, only models with very recent par- 
ticle injection (at z <J 0.01) have any significant radio emission. 
The high energy particles tend to be in steady state, and the 
radio spectra of the models tend to be quite similar. Three of 
the most divergent radio spectra are shown in Figure M The 
radio spectrum of Model 23 is typical of those for steady-state 
models with continual particle injection, and applies to most 
of the models with significant radio emission. In Model 6, the 
electrons have not yet reached steady-state, and the radio spec- 
trum is flatter. All of the electrons in Model 9 were injected at 
Zi = 0.01, and the highest energy particles have been removed 
by IC and synchrotron losses. This spectrum is quite steep. 

In general, HXR emission and radio synchrotron emission 
are expected only in clusters with very recent or current injec- 
tion of relativistic electrons. Thus, both measure the current 
rate of particle injection. For example, if the particles are ac- 
celerated in ICM shocks, HXR and radio emission would be 
expected only in clusters which are currently undergoing (or 
which very recently underwent) a merger. 

6.3. Optical and UV Emission 

One also expects that the lower energy portion of the cos- 
mic ray population in clusters will produce diffuse optical and 
UV emission. Diffuse optical emission is known to exist in 
many clusters, particularly in those with central cD galaxies 
(e.g., Boughn & Uson 1997). Although the origin is not com- 
pletely understood, the optical colors of the diffuse light sug- 
gest that it is due to old stars, which may have been stripped 
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from cluster galaxies. It is likely that the near or vacuum UV 
are better regions to detect low surface brightness diffuse emis- 
sion due to IC emission, since the older stellar population in E 
and SO galaxies (and, presumably, the intracluster stellar popu- 
lation) are fainter there. 

Table || gives the power L V (UV) at a wavelength of 2000 
A and the best-fit power-law spectral index ao P t uv between 
2000 A and 8000 A. The UV powers of most of the mod- 
els lie within a relatively narrow range of about an order of 
magnitude. Unless the electron density in the ICM is much 
higher than 0.001 cm -3 and Coulomb losses are catastrophic 
or the electron population in clusters is very old, one would 
always expect a significant population of lower energy elec- 
trons. The IC spectra in the optical and UV are fairly flat, 
with spectral indices of -0.5 ;> ceopt.uv J> 0.3. In models with 
current particle injection, the lower energy electron population 
should be approaching steady-state with the Coulomb losses 
(eq. [y]). Then, the expected power-law for the IC emission is 
a « —(p — 2)/2 w -0.15. For models with only an initial elec- 
tron population, the electron spectrum is expected to be nearly 
flat at low energies, and the spectrum is given by equation (p5[). 
Thus, the spectral index should be slightly positive in these 
models. 

The predicted fluxes are rather low when compared to the 
diffuse optical emission seen in clusters of galaxies or to the 
sensitivity of current and UV instruments for detecting diffuse 
emission. 

7. CONCLUSIONS 

Models for the integrated population of primary cosmic ray 
electrons in clusters of galaxies have been calculated. The evo- 
lution of the relativistic electrons included the effects of energy 
losses due to IC scattering, synchrotron emission, Coulomb 
losses to the ICM, and bremsstrahlung. For typical cluster pa- 
rameters, the combined time scale for these losses reaches a 
maximum of ~ 3 x 10 9 yr for electrons with a Lorentz factor 
7 ~ 300. This maximum loss time scale is comparable to the 
Hubble time or the typical age of clusters. For relativistic elec- 
trons with either much higher or lower energies, the loss time 
scale is considerably shorter than the typical age of a cluster. 

Although the models don't depend on the detailed nature of 
the source of the relativistic electrons, we assume that they are 
accelerated by shocks, either due to the formation of the cluster 
and subcluster mergers, or due to radio galaxies. We present 
models in which the electrons are all injected at some given 
time in the past (e.g., the epoch of cluster formation or the time 
when a powerful radio galaxy is active). We also consider mod- 
els with continuous electron injection, either at a constant rate 
or at a rate given by the self-similar solution for secondary in- 
fall onto clusters (Bertschinger 1985). Hybrid models with both 
initial electron populations and continual particle injection were 
also calculated. Based on the results for particle acceleration is 
supernova shocks, in most of the models we assume that the 
energy spectrum of the injected electrons is a power-law with 
an exponent of 2.3. We also consider a few models with a bro- 
ken power-law, where the break occurs at the energy at which 
protons become relativistic. 

Simple analytical solutions are given for models in which the 
electrons are either all injected at a single time in the past, or 
in which there is a constant rate of injection of electrons but 
no initial electron population. These solutions apply in the lim- 
its of either high or low energies, where IC and synchrotron or 
Coulomb losses dominate. For models in which all of the elec- 



trons were injected in the past, there is a high energy cutoff to 
the present electron distribution for 7 > 7 mflt . Expressions for 
the resulting populations and 7 mm were given which included 
synchrotron losses and the cosmological variation in IC losses 
due to the redshift of the CMB radiation field. At low ener- 
gies where Coulomb losses dominate, the electron distribution 
function M7) tends to a constant value, independent of 7. 

Analytic solutions were also given for models with constant 
electron injection. At very low or high energies, the particle 
distribution approaches steady state. If the electrons are in- 
jected with a power-law distribution, the steady state distribu- 
tion is one power steeper at high energies and one power flatter 
at low energies (e.g., Ginzburg & Syrovatskii 1964). A self- 
similar solution which connects these steady-state solutions to 
the time-dependent solution when the losses are not rapid was 
also derived. 

The numerical models for the electron population are nor- 
malized such that the total energy of relativistic electrons in- 
jected with 7 > 300 is 10 63 ergs. For most of the models, 
this leads to a present-day energy in electrons with 7 > 300 
of ^10 62 ergs. Because the equations for the electron popula- 
tion are linear, models for other values can be derived by simple 
scaling from these values. 

The most important factor affecting the energy distribution 
of cosmic ray electrons in the cluster models is the history of 
injection of the electrons. The numerical models show that only 
clusters in which there has been a substantial injection of rela- 
tivistic electrons since z <J 1 will have any significant popula- 
tion of primary cosmic ray electrons at present. For older elec- 
tron populations, the high energy electrons are removed by IC 
losses, whose rate increases with redshift. Lower energy elec- 
trons (including initially higher energy particles which have lost 
energy by IC scattering) are removed by Coulomb losses. 

In models with current (or very recent, z <C 0.1) particle in- 
jection, the electron population is large at low energies and ex- 
tends to very high energies. In models with ongoing electron 
injection, the distribution quickly reaches steady-state at high 
(7 ;> 300) and low (7 <J 300) energies. Due to IC losses, the 
high energy electron distribution is steeper than the injected 
one (by one exponent for a power-law injection spectrum). 
Coulomb losses cause the low energy distribution to be flat- 
ter than the injected spectrum (again, by one exponent for a 
power-law injection). 

In models in which there is no current or recent particle in- 
jection, the high energy electron distribution has a cutoff for 
7 > 7ma.v Thus, no very high energy primary electrons should 
be present. The low energy distribution is also considerably 
flatter than that in models with current particle injection. In 
models with a large initial population of particles, but also with 
a significant rate of current particle injection, the electron dis- 
tributions are a simple combination the behavior of the initial 
population models and the steady injection models. There is 
a steep drop in the electron population at 7„, I/V , but higher en- 
ergy electrons are present at a rate determined by the current 
rate of particle injection. The electron population at 7 < j maj! 
is fairly flat, but eventually start to increase significantly with 
decreasing energy at a lower energy. 

Increasing (decreasing) the ICM thermal gas density de- 
creases (increases) the number of low energy electrons (7 <; 
100), since these electrons are removed by Coulomb losses to 
the ambient plasma. If the magnetic field is greater than gener- 
ally expected over most of the volume of a cluster (B 3 /iG), 
synchrotron losses will reduce the number of high energy elec- 
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trons beyond the effect of IC emission. 

The portion of the electron spectrum where the different 
models exhibit the greatest level of agreement is the region near 
7 ~ 300, where the loss time of the electrons is maximum and 
where it reaches values of ~ 3 x 10 9 yr, which are comparable 
to cluster ages. The general prediction is at all clusters should 
contain a significant electron population at these energies, as 
long as there has been significant particle injection since z < 1 . 

I also calculated the IC and synchrotron emission from 
these models. In models with steady particle injection with 
a power-law exponent p, the IC spectra relax into a steady- 
state form. At low energies, the spectrum is a power-law with 
a -{p-2)/2 -0.15, while at high energies it is a power- 
law with a as — (p/2) ~ —1.15. These two power-laws meet at 
a knee at v ~ 3 x 10 16 Hz. In models with no current particle 
injection, the cutoff in the electron distribution at high energies 
(7 > Imax) results in a rapid drop in the IC spectrum at high 
frequencies. At low frequencies, the IC spectra of these mod- 
els are fairly flat or even slightly rising with frequency. The 
IC spectra of hybrid models with both an initial population of 
particles and particle injection at the present time are interme- 
diate. In models in which the current rate of particle injection 
provides a small but significant fraction of the total electron en- 
ergy, the spectra show an extended hump at low frequencies 
O £ 10 17 Hz), with a rapid fall off above v ~ 10 16 Hz. How- 
ever, they also have an extended hard tail of emission at high 
frequencies, which has a power-law spectrum with a spectral 
index of a » — (p/2) w —1.15. 

In the models, EUV and soft X-ray emission are nearly ubiq- 
uitous. This emission is produced by electrons with 7 ~ 300, 
which have the longest lifetimes in clusters. The lifetimes of 
these particles approach the likely ages of clusters, and thus 
particles at this energy are likely to be present in most clus- 
ters. This may explain why EUVE has detected EUV emis- 
sion from all of the cluster observed which lie in directions of 
sufficiently low Galactic columns (Lieu et al. 1996a,b; Mittaz, 
Lieu, & Lockman 1998; Sarazin & Lieu 1998). In the mod- 
els, the increase in the IC losses in going from the EUV to the 
soft X-ray region causes the spectral to decline rapidly in this 
region. More steeply declining spectra are expected in clusters 
without current sources of particles or with lower thermal gas 
densities. Assuming that the acceleration of relativistic elec- 
trons in clusters is due to ICM shocks, this may explain why 
the EUV spectrum in the very regular cluster A 1795 is steeper 
than that in Coma (Mittaz, Lieu, & Lockman 1998), and why 
the spectrum in A 1795 gets steeper at larger radii. 

The IC emission also extends down in frequency into the UV, 
optical, and IR bands. The spectrum in this region is expected to 
be fairly flat, with power-law spectral indices in the range from 
-0.6 to +0.3. The steeper spectra are for models with significant 
current particle acceleration, while the flatter or rising spectra 
are for models with older electron populations and/or broken 
power-law electron acceleration. In the optical, the luminosi- 
ties and expected surface brightnesses are below the level of 
diffuse stellar light observed in clusters. This emission might 
be easier to detect in the vacuum UV, where the older stars as- 
sociated with E/S0 galaxies are fainter. However, the surface 
brightnesses are rather low. The UV and EUV emission from 



clusters might make a contribution to the diffuse extragalactic 
ionizing radiation, although crude estimates suggest this will 
not be very significant. 

In much of the X-ray spectral band (from 0.5 to 10 keV), the 
dominant luminosity in most clusters is the thermal emission 
from the hot ICM, and IC emission will be difficult to detect. 
However, at energies ;>20 keV, IC emission can again become 
observable. I calculated the hard X-ray (HXR) IC emission 
from the cluster models. Radiation at these energies is pro- 
duced by high energy electrons (7 ~ 10 4 ) which have rather 
short lifetimes. As a result, this emission is only expected in 
clusters which have current (or very recent) particle accelera- 
tion. The luminosity of HXR emission is primarily determined 
by the current rate of particle acceleration in the cluster. The 
same particles which produce HXR radiation by IC scattering 
will also produce radio emission by synchrotron emission at 
observable frequencies. I have also calculated the radio spectra 
produced by the models. In general, the high energy electrons 
which produce the HXR and radio emission should be nearly 
in steady-state, and the predicted spectra for most models with 
significant HXR or radio emission are approximately power- 
laws with spectral indices of a w -1 . 1 . 

If the main source of primary relativistic electrons is accel- 
eration in ICM shocks, then HXR tails and diffuse radio halos 
should be confined to clusters which currently have strong ICM 
shocks. Thus, we would expect that these would mainly be 
clusters presently undergoing significant mergers. In the case 
of radio halos, there is some evidence that this is indeed the 
case (e.g., Tribble 1993). At the same time, large cluster radio 
halos are quite rare. This may be related to the short lifetimes of 
the high energy particles required and to the need for the merger 
to be at a stage where the shocks are very strong; alternatively, 
some other physical properties (e.g., magnetic field intensity or 
geometry, or high shock Mach number, . . . ) may be required to 
produce high energy electrons in ICM shocks. 

Because the HXR and radio halo emissions are produced by 
particles with very short lifetimes while the EUV emission is 
due to electrons whose lifetimes are a significant fraction of the 
likely ages of clusters, one would not expect the same power- 
law spectrum which fits the HXR and radio emission to extend 
to the EUV emission. In fact, EUV emission is seen from clus- 
ters with no detectable radio halo emission (such as A 1795), 
so it is clear that the EUV and HXR/radio cannot be connected 
with a single power-law spectrum with a normal spectral index. 
The EUV and HXR/radio should connect continuously in clus- 
ters where a large fraction of the cosmic ray electrons are being 
generated at the present time (i.e., in an ongoing or very recent 
strong merger). In the more general case, one expects a broad 
peak in the IC emission in the EUV band, a rapid drop-off in 
the spectrum from the EUV to the X-ray band, and HXR/radio 
emission proportional to the present rate of particle acceleration 
(Figures and ||]). 
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